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Focussing on a paradigmatic small system consisting of two coupled damped oscillators, we survey 
the role of the Levy-Ito nature of the noise in the thermal conductance. For white noises, we 
prove that the Levy-Ito composition (Lebesgue measure) of the noise is irrelevant for the thermal 
conductance of a non-equilibrium linearly coupled chain, which signals the independence between 
mechanical and thermodynamical properties. On the other hand, for the non-linearly coupled case, 
the two types of properties mix and the explicit definition of the noise plays a central role. 
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I. INTRODUCTION 

The law of heat conduction, or Fourier's law, i.e., the 
property by which the heat flux density is equal to the 
product of the thermal conductivity by the negative tem- 
perature gradient pj is a paradigmatic manifestation of 
the ubiquitous laws of thermodynamics 0- Recently, it 
has stoked a significant amount of work on its explicit 
derivation for large Hamiltonian systems 0-H|. In this 
context, models with anharmonic coupling succeed in dif- 
fusing energy, but the analytic solutions thereto are very 
demanding, even for the few cases where that is possible. 
Since they allow a larger number of exactly solvable cases, 
small systems are worthwhile [f| and particularly relevant 
in chemical physics and nanosystems Q ■ In the scope of 
analytical methods, we highlight the time averaging of 
observables that endure a stationary state [H, Q . This 
account has several advantages, namely compared with 
the Fokkcr-Planck approach, which cannot be applied to 
those cases where higher than second order cumulants 
of the noise are non-vanishing and also significant. This 
comprises Poissonian [Iol[l3lll4| and other non-Gaussian 
massive particles [ll| as well as other cases where the in- 
teraction with a reservoir is described by a process with 
a non-zero singular part of the measure when a Levy-Ito 
(LI) decomposition is applied fl2j . 

Stemming from these facts, we perform a time averag- 
ing study of a small non-equilibrium system composed of 
two damped coupled oscillators at distinct temperatures 
and determine the explicit formula of the Fourier's law 
for linear and non-linear cases. In spite of its simplicity, 
the former has relevant traits: i) it is a non-equilibrium 
system; ii) Its heat flux definition is well known; in) It 
is adjustable to different kinds of reservoirs; iv) It can be 
expanded into a infinite chain with a nearly direct appli- 



cation of the results of a N = 2 block; v) It represents 
the result of Langevin coloured noises by a renormalisa- 
tion of the masses [§[ and vi ) Linearity is still a source of 
important results in many areas [15l4lq . 



II. MODEL 



Our problem focus on solving the set of equations, 



dvj (t) 
1 alt 



= -kxi (t)-7«j (*)~2J fc 2Z-i [xi (t) - xj (t)] 
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1 = 1 



(1) 

with Vi (t) = , where <E {1,2} and k\ and k% 

are the linear and non-linear coupling constants, respec- 



tively. The system is decoupled (linear) for k 



1(3) 



0. 



The transfer flux, j 12(f), between the two particles reads, 

Jl2 (t) EE - J2 ^ [*1 (*) - ^ (*)] 2 ' _1 [«1 (*) + ^2 (t)\ • 
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The term, r\i (t), represents a general uncorrclated 
Levy class stochastic process with cumulants, 

(77^ (ti) ... r) in (t n )) c = A(ti,n) S n i2 ... S ln _ ± in x 

5(ti-t 2 )...5(t n - 1 -t n ).(3) 

From [l9j . we either have two or infinite non-zero cu- 
mulants. The former corresponds to the case in which 
the measure is absolutely continuous, characterising a 
Brownian process. In Eq. ©, A(t, n) is described by the 
noise; If it is Wiener-like, W(t) = r](t')dt', A(t,n) 
is time-independent and equal to a 1 for n = 2 and zero 
otherwise (a is the standard deviation of the Gaussian). 
Among infinite non-zero cumulant noises, we can include 
the Poisson process for which A(t, n) equals $™ X(t) [2fj| . 
with <!> being the p(&) independent and identically dis- 
tributed magnitude and X(t) the rate of shots. Herein, 
A is time-independent without loss of generality. For 



-Vi (t) 



2 



ki = k^ = 0, Eq. (JTJ) is totally decoupled and the so- 
lutions to the problem of homogeneous and sinusoidal 
heterogeneous Poisson noises can be found in Ref . . 



III. RESULTS 



Laplace transforming Xi(t) and Vi(t) we obtain, 



®i ( s ) — R(s)^3 ( S ) + 



^ylima-M)/// 

SXi(s) = Vi(s) 



n 



[Zi (i q n +a) — Xj (i g n +a)] 



s-(igi+i(j 2 +ii?3+3a) 



(4) 



r 



(Re(s) > 0) with H(s) = (ms 2 + 7 s + fc + /ci). The 
solutions to Eq. (0} arc obtained considering the rel- 
ative position, vb (s) = x\ (s) — .T2 (s), the mid-point 
position, fs (s) = (5i (s) + X2 (s)) /2, as well as the re- 
spective noises f\D (s) = ryi (s) — 772 (s) and 7)5 (s) = 
(771 (s) + 772 (s)) /2. After some algebra it yields, 



with i?'(s) 



#H lim f + °° 11 
fl'(s) Q _).o J -° 
s = Vs(s) 
' S R"(s) 



= 1 4!f fpiiqi+a) 



(ms 2 +js + k + 2ki) and R"(s) 



(5) 



(m s 2 + 7 s + fc) . Reverting Eq. [5] we get X\ (s) and 



5^2 (s). Concomitantly, wc must compute the Laplace 
transforms of 771 and 772, 



fail 0i) ■■• *7i„ (2«)> c = / °° ft dt h exp -E"=i^^ 

X fail (*l) ■ • • Vin (*n)) c 



that arc employed in the averages over time 



^ii Z2 ■ ■ • ^i n _i i n • 



(6) 



z->0 



= lim 

Z,£->0 



=limz/ / / 6(t-t!) 6(t-t 2 ) e~ zt (x™(h) vZ{t 2 )) c dt 1 dt 2 dt 



00 m+n 



n dqi 
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n? = i(ig^+g) (n£T£(i« + e 
z - CrST i ft' + ( m + 



(7) 

(8) 



Allowing for a contour that goes along the straight line 
from —p + i e to p + i e and then counterclockwise along 
a semicircle centered at + ie from p + ie to — p + ie 
(p — > 00 and e — > 0), we realise that two situations oc- 
cur: cither the calculation of the residues leads us to a 
term like —2— u. with (u. w) =fi 0, which vanishes in the 
limit z — ► 0, or to -u, which is non-zero. The prob- 
lem solving now resumes to an expansion of Eq. ([2]) in 
powers of k$, actually the expansion is carried out in 
powers of k^T/k 2 [2l[ (see Appendices). Hereinafter, 
besides the results of the Brownian thermostats, which 
corresponds to the pure continuous measure, we make 
explicit the conductance for the Poisson reservoirs, the 
epitome of singular measures. Nevertheless, the results 
for other non-Gaussian noises can be obtained following 
our methodology yielding the same qualitative results. In 



first order the transfer flux reads, 



with, 



and 



7 W 
J 12 



As) 
J12 



W) = (i2 } ) + (ig } ) + + O (k 2 ) , (9) 



• (0)\ _ kj [A 1 (2)~A 2 (2)] 
J 12 / 4 mfej+7 2 (fe+fci) 



3 t , (2k+k 1 )[A 1 (2Y~A 2 {2Y 

a 7 K l K 3 



(fe+2fei) 7 2 (fe+fci)+m fc 2 



7 ^W (W2)J _^ (2)J]) 
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A 2? 



(10) 



(11) 



for the Poisson case and ^Jx2 ^ = f° r the Brownian case 
(see Appendices). For Poisson, when A -1 < 1 (keeping 



3 



the temperature fixed), the weight of the singularity of 



the noise measure dwindles and \ii2 / — * 0- The coeffi- 
cients in Eq. (fTTJ) are, 

7V = 7 2 (5 fc + 3/ci)+m(3fci+4fc 2 + llfcfci) , (12) 
V= [m(4fc + 9fci) 2 + 67 2 (2fc + 3/si)j x (13) 
[3 7 4 + m 2 fc 2 +4m 7 2 (fc + A: 1 )] [ 7 2 (fc + fci)] . 

Thence, we are finally in the position to compute the 
thermal conductance, 
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dAT 



AT ' 



(14) 



Ty-T 2 

Resorting to single particle results and the equipartition 
theorem [1(| , we relate the cumulants of the noise and the 
proper temperature, Tj, namely, Ai(2) = 2jT i: yielding 
a thermal conductance, k = + + + O (kf). 
Equations l[9"])-(|lip pave the way to the following as- 
sertion; When interacting particles are subject to white 
reservoirs and coupled in a linear form, the explicit ther- 
mal conductance is independent of the specific nature of 
the noise, namely the outcome of their Levy-Ito decompo- 
sition, whereas for non-linear coupling the nature of the 
measure of the noise (its decomposition) is pivotal. In 
other words, the linear case is heedless of the measure of 
the reservoirs and it only takes into consideration their 
temperatures. Hence, is exactly the same, either we 
have a Wiener noise (continuous measure), which is the 
standard noise in fundamental statistical mechanics stud- 
ies 0, [22| , or a Poisson noise (paradigmatic case of sin- 
gular measure). Although the linear coupling result has 
only been explicitly proved for two particles, it is valid 
for general N. In fact, for a linear chain, the local energy 
flow can be written as a function of the cumulants 

(Vi(z)m(z')) c = ZjTi and (rn +1 (z)r]i +1 (z')) c = 2jT. l+1 , 
wherein the dependence on the specific nature of the noise 
is eliminated, except for the respective temperatures. On 
the other hand, if the nature of the noise affects the con- 
ductance of the simplest coupling element, the conduc- 
tance for generic chains is also changed. This result is un- 
expected since contrarily to single particle linear cases, 
wherein the LI nature is already relevant [HI El, for 
coupled systems the LI composition does solely become 
significant when the interaction between the elements of 
the system happens in a non-linear way. Only in this 
case higher-order cumulants of the noise, which can be 
learnt as higher-order sources of energy, influence the re- 
sult. For the same (fci,^), in decreasing the singularity 
by soaring A, the two thermal conductances tally. 

To further illustrate these results, we have simulated 
cases of equally massive particles subject to Wiener and 
Poisson noises at different temperatures, T. For the 
former, we have T = <r 2 /2 whereas for the latter, we 
have assumed a homogeneous Poisson process with a 



rate of events A, with a random amplitude, $, expo- 
nentially distributed, p{<&) ~ exp[— $/$], which yields 
T = A < I ,2 /7 [Tol| - In Fig. [TJ we depict linear coupling. 
It is visible that after a transient time, t* , the system 
reaches a stationary state and (712) becomes equal to 
(712)7 whatever the reservoirs. In fact, even more com- 
plex models, such as linear chains of oscillators, verify the 
k = kS®> property. Still, this is valid when each particle 
is perturbed by different types of noise, e.g., a Brownian 
particle coupled with a Poissonian particle. The instance 
where the noises are of different nature gives rise to an 
apparent larger value of the standard deviation^! 
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FIG. 1. (Colour on-line) Average exchange flux (j'12) of a two 
massive particle system for different combinations of paradig- 
matic types of noise with T\ = 10 and T2 — 121/10, m = 10, 
7 = k — 1, fci = 1/5, &3 = and A = 10 for Poissonian 
particles. After the trasient, k agrees with the theoretical 
value, k = 21/800 = 0.02625, with the fitting curves lying 
within line thickness. The averages have been obtained aver- 
aging over 850 x (5 x 10 5 ) points. The discretisation used is 
St — 10~ 5 with snapshots at every At = 10~ s . 



In turning ks 7^ 0, the composition of the measure of 
the reservoirs comes into play. In Fig. [21 we show the dif- 
ference between equivalent Brownian and Poissonian par- 
ticles with a good agreement between the averages over 
numerical realisations and the respective (first order) ap- 
proximation. For the same temperature, the larger the 
value of k, the larger the value of defining a 10% dif- 
ference between numerical values the approximation. 



IV. FINAL REMARKS 

To summarise, we have studied the thermal conduction 
in a paradigmatic mechanical system composed of two 
coupled damped harmonic oscillators subject to generic 



1 Although the computation of &^ 12 = (^12) — O12) 2 is possible, 
we have set it by as it demands a mathematical tour de force 
likely to yield a lengthy formula with little grasping information. 
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FIG. 2. (Colour on-line) Comparison between numerically 
obtained values (symbols) and the first order approxima- 
tion of thermal conductance from Eqs. (|9))- (|11[) for different 
temperatures pairs, namely A = {10,^}, B = {10, 
C = {10, 2^} with m = 10, 7 = k = 1, ki = 1/5 and A = 1 
for Poissonian particles. 



noises, which can be understood as a concise way to de- 
scribe non-equilibrium problems. By averaging in the 
Laplace space, we have been able to determine the con- 
ductance of a linearly coupled system and approximate 
formulae for non-lincarly coupled particles. We have 
shown the conductance of the former is independent of 
the nature of the (white) noise, namely its Levy-Ito de- 
composition structure. This result is unexpected since 
the measure of the thermal bath plays a major role for 
single particle properties. The dependence on the noise 
only emerges when there exists transfer of energy in a 
non-linear way and higher-order cumulants of the noise 
enter in the calculations. In the case of Poisson noises, 
we show that the difference to Brownian noises becomes 
negligeble when the ratio between the coupling constants 
and the rate of events is small. Our calculations evidence 



the independence of the thermodynamical properties of 
the system from the nature of the reservoirs in linearly 
coupled systems. On the other hand, when the coupling 
is non-linear, the nature of the reservoirs affects the con- 
ductance, which represents a mixture between mechanic 
and thermodynamical properties of the system. 

Our results have direct implication for the study of 
the thermal conductance of systems under the influence 
of noises other than Wiener, for instance: i) solid state 
problems wherein shot (singular measure) noise is re- 
lated to the quantisation of the charge [2j|; ii) RLC 
circuits with injection of power at some rate resem- 
bling heat pumps [lH; in) Surface diffusion and low 
vibrational motion with adsorbates, e.g., Na/Cu(001) 
compounds [25j; iv) Biological motors in which shot 
noise mimics the nonequilibrium stochastic hydrolysis of 
adenosine triphosphate [l4[ and v) Molecular dynam- 
ics when the Andersen thermostat is applied. Actu- 
ally, in molecular dynamics [24j . the Langevin reservoir 
is just one in a large collection of baths represented by 
our definition of noise ((3]). In these problems, for non- 
lincarly coupled elements the experimentally measured 
energy flux will be greater than the energy flux given by 
Langevin reservoirs at the same temperatures and equal 
if coupling is linear. At the theoretical level, the method 
is worth being used to shed light on non-linear chains as 
well. Within this context, the feasible approach is once 
again to consider a perturbative expansion of the non- 
linearities in the problem. 
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Appendix A: Heat flux definition 

As indicated in the main document the model consists of two masses connected by linear and non-linear springs so 
that, 



x\ = Vi, 

X 2 = V 2 , 

mil = — kxi — ki(x± — x%) — k^{x\ — x 2 ) 3 — 7±i + 7?i 

mx 2 = -kx 2 - h(x 2 - x\) - k 3 (x 2 - X1) 3 - -yx 2 + f] 2 



3 (Al) 



The flux of energy J±^ 2 defines the thermal conductance. For this unidimensional mechanical system, the trans- 
mitted power is by the the instantaneous power difference 0]: 

dWi^ 2 - dW 2 ^i Fi^ 2 v 2 - F 2 ^ivi 
Jl ^ 2 2dt 2 ' (A2) 
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where F a ^p is the force exerted on particle (3 by particle a. 

Assuming a conservative potential, the form of the force between the particles will be, 

F a _>j3 = -ki(x a - xp) - k 3 (x a - xp) 3 . (A3) 

Our goal is to derive a systematic expansion for the time-average ji 2 = (Ji^ 2 ) = K (T 2 ~T\), where k is the thermal 
conductance of the model. Thus, we can write 

ii 2 = ^ fi^jWa ) 

(A4) 

where the terms of the form < x™vi > and < x 2 l v 2 > vanish identically as t — > 00 0. 

1. The Gaussian case 

The noise functions are assumed to be white and Gaussian, 



< r/i(ti) > c = 0, 
<m(ti)m(ta)>c = 2jT 1 S(h-t 2 ), 

< m{ti) >c = 0, 

<m(h)m(t2)>c = 2 1 T 2 5{t l -t 2 ), 



«i(s) = sii(s), 
v 2 (s) = sx 2 (s), 



using i?(s) = (m s 2 + 7 s + (fc + fci)) = ?n(s 2 + 0s + cj 2 ) = (s — C+)( s — C-)j where 



we have, 



and 



*i« = ife saw + is - ifo n-^o rr & /_ ~ & /_r & >< 

(xi(igi+q)— a2(»gi+a)) 0ri(»g2+t»)— £2(»g2+oQ) (xi(iq3+a)— Z2(»g3+o0) 
s— (iiji +i<j2 +i<?3 +3a) 

saw = t^j + is + t^i Bmo^o rr & & /r 3? >< 

(£i(igi+q)— a 2 (»gi+a)) (£1 (iq 2 +oQ— £2 (»g2+oQ) (xi(iqs,+a)-x 2 (iq3+a)) 
s— (iqi +iqz+iq3 +3a) 

An straightforward way to write a recurrence equation for the problem above is to take the difference: 

(1 + ife) tow - **(■)) - - s& ii-^o rr & /r /_» if >< 

— £ 2 (igi+a)) (xi (iq2+oQ — £2 (i<?2+oQ) (£1 (»g3+oQ -£2 (ig3+a)) 
s — (igi +i<?2 +iq~3 +3a) 

Defining i?'(s) = (ms 2 + 7 s + (fc + 2fci)) = (s 2 + 6»s + ujj) = (s - Ci+)( s - Ci-), where, 



(A5) 



The initial conditions will be assumed to be 

a;i(0) = «i(0) = 13(0) = wa(0) = 0. (A6) 

The Laplace transforms read, 



(A7) 



C±--2±|V4o; 2 -0 2 (A8) 



(A9) 



(A10) 



(All) 
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we obtain the difference, 

r D = xi(s) - x 2 (s) 

m(s)-il2(s) 2fc 3 i- r+oo dqi r+oo d q 2 r+oo d q 3 , A1 „, 

(x 1 (iq 1 +a)-x 2 (iqi+a)) (ii (192 +a) -x 2 (iq2+a)) {Zl (»93 +a) -x 2 (iq-j+a)) 
s-(i qi +iq 2 +iq 3 +3a) 

Defining R"(s) = (m s 2 + 7 s + k) = (s 2 + 6»s + w|) = (s - C2+)(s - C2-), where, 



we obtain the sum as, 



Inverting the relations it yields, 



~ _ x 1 {s)+x 2 {s) 
' b 2 



2 R"{s) ■ 



xi(s) = rs + ^f-, 



X2(s) = 

In the same way, we can define the difference and the average of the noise as, 

Vs(s) = fcife 

fj D (s) = m - m- 

We can now express the recurrence relations for the new variables, 

_ VP (s) _ 2fc 3 i- r+°o d qi r+oo dq 2 r + oo dq 



x 

00 Itt 



f D (iqi+a) f D (iq 2 +a) fiAiqz+a) 
s—{iqi +192+193+30) 



The Laplace transforms of the noise are, 

< VdMvdM >c=- 2 ^ (Tl+T2) 



The series expansion for the thermal flux reads, 

2 



(A13) 



(A14) 



(A15) 



(A16) 



< ^(«i)%(* 2 ) > c = (A17) 



Appendix B: Heat conductance 
1. General expression 



312 = ~ (Oi - x 2 )(v 2 + Vl )} - -| ((X! - x 2 f(v 2 + (Bl) 



We can write (a, b = 1,2), 

"6 ) c = lim^o lim e ^o JILi ^ .Co IlILi if z -(^r =1 «ft+(m+n) 6 ) (nr=i ^(^fc + e) FGU + e )) c 



lim^o nr=r ^ z _ (Er + r 4 +(m+n)£) n?=i(v» + <o (nr=i n + <o) c • 

(B2) 
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We can thus express the heat flux (j'12) as, 

(J12} - -ki lim 2 ^o lim e ^o JZo I-oo z -( igi + iq2+ 2e) (*& + e ) + e ) fs(*fc + e)) c - 

x {(iq 4 + e) (r D {iqi + e)r D (iq2 + e)f D (iq 3 + e)rg(ig 4 + e)) c } . 
We now proceed to expand the heat flow in powers of k 3 . 



2. Order zero on £3 



The order zero term is, 



J-. v- _ _, roo den f 00 dq 2 z{iq 2 +e) (fjp (jgi +e) fjs (»92+e)) c 

U12/o — (.imil^mile^o J_ m 2lr J.^ 2 tt z-(igi+j 92 +2e) i?/(i 9 i+e) i?"(ig 2 +e) 



-00 2tt J- 00 27r z — {iq\+iq 2 +2e) R' {iqi +e) i?" (?'q 2 +e) 

7(ri-T a )fci i- i- r°o poo dg 2 z (»g2+^) 1 

m 2 w^-z-tO llui e->-0 J_ oc 27r J_ oc 27r 2 -(,:<j 1 +ig 2 +2e) j?/(i<ji+e) R"(iq 2 +e) {iq 1 +iq 2 +2c) 



which is already a well known result 0, ■ 



1 7 (T 1 -r 2 )fc 2 

2 7 2 (fc+fc 1 ) + fc 2 m ' 



3. Order 1 on &3 



There are two contributions to order 1 on k 3 : one from the quadratic term and another from the quartic one. The 
quadratic term reads, 

(j^ 12 = kl lim^o Hnie^O lim Q ^ f °° ... f + °° ^ ^gHo , B „ 2 \ , <^(W3+a)f D (ig 4+ a)f J (i g!i +a)f s (i ga+ 0> r 



i<?l+c-(i<?3+*94+i<?5+3a) 



12 7 2 fc! fc 3 lim 2 ^ lim^o hm^o ^ ■ • • /_ 



oo r+oo 2 (tg 2 +£) 



X757 



00 27r ' J — 00 2iz z — {iq\+iq 2 +2e) iqi+e—(iq^+iq4-\-iq^+3a) 

(Ti+T 2 ) (T1-T2) 



X 



fl'(igi+e) ii'(iq3+a) #'(194+01) -R' («<?5 +a) -R"('i<?2 + e) (i<?3+i<?4+2o) (('(/.-. +<i +iq 2 +f-) 



3 (^fc-mfc^Tx'-T., 2 ) 
2 ' 3 (fe+2fei)[mfc 2 +7 2 (fe+fci)] 



The quartic term reads, 

ij^; = -k 3 lim^o iim e ^o n, • ■ • r — 



(B4) 



00 27r " ' 'J— 00 27r 2 — (2gi+iq2+^3+iij4+2e) 

x {(ig 2 + e) (f D {iqi + e)fs(iq2 + e)f D (iq3 + e) + e))J 



-6 7 2 A;3 lim^o lim e ^o f ■ • ■ f °° r I? I + ^ * ran 

I o z—ru t— fu J^oo 27r J-00 27T z~(iq 1 +iq 2 +iq 3 +iq i +2e) (B5) 



x (jg2+e) (T1+T2) (Ti-T 2 ) 

fl'(igi+e) /?'(ig 3 +Q) -R'(igi+c) -R"(ig 2 +£) (iq;- i +iq i +2a){iq 1 +iq 2 +2t) 

_ 3 Tfci fc 3 (ri 2 -r 2 2 ) 

2 (fc+2fci) [7 2 (fc+fci)+mfc 2 ] ' 



The total term reads, 



/ 2 (fc + 2/ci)[m/c 2 +7 2 (fc + fci)] 2 
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By the relation j\i = —kAT, we obtain the expansion of the thermal conductance as, 

(B7) 



iK , 3 7 J fci fc 3 (2fc+fc 1 )(T 1 +T 2 ) n(T,2\ 

7 2 (fc+fci) + fc 2 m ~r 2 (fc+2fci)[mfc 2 + 7 2 (fc+fci)] 2 " t " ^^3^ 



Appendix C: Poisson Conductance 

An interesting use of our formalism is to compare the result of Eq. (|B7[) with the case of Poisson noise. For Poisson 
we have [Ioj |. 

< V(S1)V(S2) >c = 

(CI) 

< 77(si)7?(s 2 )^( S3 )f?(s4) > c = { Sl l 4 S2 \* +Si , 

where the expression for the fourth order cumulant can be checked by dimensional analysis and compared with that 
of reference [1(| • 

The combinations < fjs,Dfjs,D > need to be reexamined. For the Poisson noise, they yield, 

< fj D ( Sl )a D (s 2 ) > c = 27 f;: 2 r2) , 

< Usi)vsM >c = 

(C2) 

<vsMvdM>c lS §^r> 

<Vs(si)fj D (s 2 )v D (s 3 )fj D (s 4 )> c = i^gS^. 

The zero-th order on fc 3 is exactly the same as the Gaussian case. The first order can be illustrative and we shall 
calculate it in the following. The quadratic term reads, 

Ta \ U 11™ U™ U m r°° d H r°° d <X> * (jg2+e) 2 fc 3 (fp(iqi+a) r D (ig 4 +a) f D (iq 5 +a) f s (iq 2 +e)) c 



-oo oo 27r z— (iqi+i</2+2e) R'(igi+e) iqi+e — (293+294+295+30) 
oo dqi roo dq 5 z (»?2+e) 



24 A 7 2 fci £; 3 lim z ^ lim e ^ lim a _> J!!^ ^ • • ■ /_ 



OO oo 2-7T 2 — (iqi+iq2+2e) 291+e — (293+294+^95+30) 



X 



x (ift + e) j-. 



(Ti+T 2 )(Ti-T 2 ) 
93 +i<J4 +2a+i(ji +i(j 2 +2e) 



2 fcife 3 (8fe 2 7 2 m+207 2 fcmfei+6fe7 4 -97 2 mfei-4fcm 2 fei-9m 2 fcJ)(Ti 2 -T 2 2 ) 
' A [7 2 (fc+fei)+mfe 2 ][6 7 2 (2fc+3fci) + 16fc 2 m+72Tnfefei+81mfe 2 J(3 7 4 +4 7 2 fcm+4 7 2 mfci+m 2 fc 2 ) 

(C3) 

The quartic term reads, 

0'ia>i4 = - k 3 lim ^o lim ^o S-co W • • • I-oo ^ z-{ lqi +lql%l+ lqi +2,) ( f 'oim + e) rs(*«2 + e) ^(^3 + e) ^(^4 + e)) £ 
= - 12 A 7 2 fc 3 lim^o lim^o ^ • ■ • I- ^ 



X 



oo 2t z-(iq% +iq2 +193+294 +2e) fl'(iqi+e) -R' (i(j3+ct) -R'('i(?4+e) -R"(ig 2 +e) 
(Ti+r 2 )(Ti-T 2 ) 



R'{iqi+e) R' (iqz+u) R' (iq&+a) R'(iqs+a) -R"(ig 2 +e) (293+194+295 +3o+29 2 +€) 



2 fci fe 3 (8fc 2 7 2 m+207 2 femfci+6fc7 4 -97 2 mfc 2 -4fcm 2 fc 2 -9m 2 fcJ)(Ti 2 -T 2 2 ) 
' ^ [7 2 (fe+fci)+?nfe 2 ][6 7 2 (2fc+3fci) + 16fe 2 m+72mfcfci+81mfc 2 ](3 7 4 +4 7 2 fcm+4 7 2 mfci+m 2 fe 2 ) ' 

(C4) 

The total Poisson conductivity at order 1 on fc 3 reads 

1 fk 2 I 3 7 3 fci fc 3 (2fc+fci)(Ti+T 2 ) 

KPmsson ~ 2 7 2 (fc ! + fe)+fe 2 m + 2 (fc+2 fc 1 )[mfe 2 + 7 2 (fc+fe 1 )] 2 

+54 7 4 ^ $ (Ti + T 2 ) + O (fcf) 
K (o) + K (D + + o (Jfcf ) 
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where 

AT = 5k^ 2 + 3j 2 k 1 +3mk 2 + 4k 2 m + llmkk 1 

V = [ 7 2 (fc + k l )+mk 2 } [67 2 (2fc + 3fci) + 16/s 2 TO + 72mfcfci +81m£; 2 ] x ^ 

x [37 4 + 47 2 m (k + fci) + m 2 k 2 ] . 
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